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We present the results of a long-look monitoring of 3C 273 with RXTE between 
1996 and 2000. A total of 230 observations amounts to a net exposure of 845 ksec, 
' with this spectral and variability analysis of 3C 273 covering the longest observation 

, period available at hard X-ray energies. Flux variations by a factor of 4 have been 

detected over 4 years, whereas less than 30% flux variations have been observed for 
individual flares on time-scales of ~3 days. Two temporal methods, the power spec- 
trum density (PSD) and the structure function (SF), have been used to study the 
variability characteristics of 3C 273. The hard X-ray photon spectra generally show 
Jjj | a power-law shape with a differential photon index of T~ 1.6±0.1. In 10 of 261 data 

segments, exceptions to power-law behaviour have been found: (i) an additional soft 
excess below 4 keV, and (ii) a broad Fe fluorescent line feature with EW ~ 100—200 
eV. Our new observations of these previously reported X-ray features may imply that 
3C 273 is a unique object whose hard X-ray emission occasionally contains a compo- 
K% ' nent which is not related to a beamed emission (Seyfert like), but most hard X-rays 

are likely to originate in inverse Compton radiation from the relativistic jet (blazar 
like). Multi- frequency spectra from radio to 7-ray are presented in addition to our 
RXTE results. The X-ray time variability and spectral evolution are discussed in the 
framework of beamed, synchrotron self-Compton picture. We consider the "power bal- 
ance" (both radiative and kinetic) between the accretion disk, sub-pc-scale jet, and 
the 10 kpc-scale jet. 

Key words: quasars: individual (3C 273) — radiation mechanisms: nonthermal — 
X-rays: galaxies 



1 INTRODUCTION 

As the brightest and nearest [z = 0.158) quasar, 3C 273 
is the ideal laboratory for studying active galactic nuclei 
(AGN). Studies of this source are relevant to all AGN 
physics, as 3C 273 displays significant flux variations, has a 
well- measured wide-band spectral energy distribution, and 
has a relativistic jet originating in its central core (see Cour- 
voisier 1998 for a review). VLBI radio observations of the 
pc-scale jet have shown a number of jet components mov- 
ing away from the core at velocities apparently faster than 
the speed of the light (e.g., Pearson et al. 1981; Vermeulen 
& Cohen 1994). The collimated jet structure extends up to 



~ 50 kpc from the core. Since 3C 273 is bright at all wave- 
lengths and on various scale sizes, it provides a valuable 
opportunity to probe the most inner part of the accretion 
disk (~ 10~ 4 pc) as well as the large scale jet (~ 10 4 pc) at 
the same time. 

3C 273 is generally classified as a blazar and is also a 
prominent 7-ray source. It was the only extra-galactic source 
of gamma-rays identified in COS-B observations (Swanen- 
burg et al. 1978), and was subsequently detected at energies 
in the 0.05 MeV to lOMeV range with OSSE (McNaron- 
Brown et al. 1995), in the 0.75 MeV to 30 MeV range with 
COMPTEL (Schonfelder et al. 2000), and above 100 MeV 
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in numerous EGRET observations (Hartman et al. 1999). 
EGRET observations helped establish that the overall spec- 
tra of blazars (plotted as vF v ) have two pronounced con- 
tinuum components: one peaking between IR and X-rays, 
and the other in the 7-ray regime (e.g., Mukherjee et al. 
1997). The low energy component is believed to be pro- 
duced by synchrotron radiation from relativistic electrons 
in magnetic fields, while inverse-Compton scattering by the 
same electrons is thought to be dominant process respon- 
sible for the high energy 7-ray emission (Ulrich, Maraschi, 
& Urry 1997). The radiation is emitted from a relativistic 
jet, directed close to our line of sight (e.g., Urry & Padovani 
1995). 

3C 273 is no exception to this picture, but an additional 
Big Blue Bump (hereafter BBB ) dominates the optical- 
soft-X-ray emission (see Paltani, Courvoisier, & Walter 1998 
; Robson 1996 for a review). Interestingly, although simi- 
lar excesses have been observed in the optical-UV region of 
Seyfert galaxies, they have not been reported for any blazar 
other than 3C 273. Although its origin is still far from being 
understood, it has been proposed that BBB may be due to 
thermal emission from the surface of a standard accretion 
disk (Shields 1978), including optically thin parts of the disk 
and a corona (see Czerny 1994 for a review). Various models 
have been suggested (e.g., Courvoisier & Clavel 1991), but 
it is generally agreed that the BBB is nearly isotropic emis- 
sion from the vicinity of the central black hole, presumably 
from the accretion disk. 

The presence of a fluorescent emission line at 6.4 keV 
is a signature of X-ray reprocessing by cold material. There 
is some evidence for the presence of a weak line at this en- 
ergy in the X-ray spectrum of 3C 273. One of the Ginga 
observations showed evidence for the line at 99% level, but 
the other observations provided only upper limits (Turner et 
al. 1990). A line detection is also reported by Grandi et al. 
(1997) and Haardt et al. (1998) in a BeppoSAX observation 
in 1996. ASCA provided only upper limits in a 1993 obser- 
vation (Yaqoob et al. 1994), whereas a broad line feature 
was clearly detected in 1996 observations. This somewhat 
confused situation probably means that both thermal (as 
for Seyferts) and non-thermal (as for blazars) emission pro- 
cesses are taking place in this particular object. However, 
it is completely unknown (i) which process dominates the 
radiation and (ii) how often the Fe line is "visible" in the 
photon spectrum. 

Many 7-ray blazars, including 3C 273, have shown large 
flux variations either, on time-scales as short as a day for 
some objects (e.g., von Montigny et al. 1997; Mukherjee et 
al. 1997). However, there are only a few blazars whose vari- 
ability characteristics are well studied in the high energy 
bands (both X-ray and 7-ray bands). For example, recent 
X-ray studies of Mrk 421, the proto- typical "TeV emitting" 
blazar, have revealed that (i) the variability time scale is ~ 
1 day, and (ii) the flux variation in the X-ray and the TeV 7- 
ray bands is well correlated on time-scales of a day to years 
(see Takahashi et al. 2000; Kataoka et al. 2001). Clearly, 
these observations provide important clues to understand- 
ing jet physics, and potentially for discriminating between 
various emission models for blazars. Although 3C 273 is a 
bright object and particularly well-sampled, it has not yet 
been possible to undertake such studies since data from pre- 



vious X-ray satellites were too sparse (e.g., the 13 observa- 
tions spanning over 5 years of Turner et al. 1990). 

In this paper, we analyze the archival hard X-ray data 
obtained with RXTE between 1996 and 2000, with a to- 
tal exposure of 845 ksec. This report of both the temporal 
and spectral variability of 3C 273 is thus based on the high- 
est quality and most densely sampled data in this energy 
band. The observation and data reduction are described in 
§2. Two temporal methods are introduced; the power spec- 
trum density (§3.1) and the structure function (§3.2). The 
X-ray spectral evolutions are summarized in §§4.1 and 4.2. 
Multi-frequency spectra are presented in §4.3. In §5, we dis- 
cuss scenarios which systematically account for the hard X- 
ray variability and spectral evolution of 3C 273. Throughout 
the paper, we discuss the energetics between the central en- 
gine and the relativistic jet, taking into account observations 
of the large scale jet by Einstein, ROSAT and Chandra. 
Finally, in §6, we present our conclusions. 

Throughout this paper, we adopt Ho = 
75 kms -1 Mpc -1 and go = 0.5. The luminosity distance to 
the source is c?l = 2.02 x 10 27 cm. 



2 OBSERVATION AND DATA REDUCTION 

3C 273 was observed 230 times with the X-ray satellite 
RXTE between 1996 February and 2000 February, with 
a net exposure of 845 ksec. The observations are summa- 
rized in Table 1. All RXTE observations were performed 
with Good Xe-16s plus Standard 1/2 modes for the Propor- 
tional Counter Array (PCA; Jahoda et al. 1996). The source 
counts were extracted from three Proportional Counter 
Units (PCUO/1/2) for the 1996-1999 observations, which 
had much larger and less interrupted exposures than those 
for PCU3 and PCU4. After 1999 May, PCU1 experienced 
the same "discharge problems" as for PCU 3 and 4, and was 
often operated with a reduced high voltage supply. We thus 
use the data from two PCUs (PCUO/2) for the analysis of 
1999 May to 2000 Feb. data (119 ksec of P40176). 

We used only signals from the top layer (X1L and 
X1R) in order to obtain the best signal-to-noise ratio. 
Standard screening procedures were performed on the 
data, using analysis software package HEASOFT 5.0 pro- 
vided by NASA/GSFC. Backgrounds were estimated using 
pcabackest (version 2.1b) for the PCA, and subtracted from 
the data. We have not used data from the High Energy X- 
ray Timing Experiment (HEXTE) on board RXTE, for two 
reasons; (i) the typical exposure for RXTE observations was 
too short to yield meaningful hard X-ray data of this variable 
source above 20 keV, and (ii) calibration problems make the 
analysis results quite uncertain. Similar problems exist for 
the All Sky Monitor (ASM) data on board RXTE. Given 
the large systematic errors already known (e.g., enhanced 
noise at certain solar angles) in ASM faint source data, we 
cannot justify the use of ASM data for detailed temporal 
studies. 

In order to obtain the maximum photon statistics and 
the best signal-to-noise ratio, we selected the energy range 
2.5—20 keV for the PCA. The overall light curves are shown 
in Figure 1. The ASM light curve (2-10 keV rates), in 10-day 
bins, is shown for comparison with the PCA flux variations 
(though one finds no clear correlation between them, be- 
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Table 1. Observation log of 3C 273. 



Obs. timet a Exposure (ksec) Data 



Obs.ID^ Obs.datc 

P1O33O0 1996 Jul 16 23:15-1996 Jul 18 13:21 

P10354 1996 Feb 02 01:12-1996 Aug 24 01:30 

P20349 1996 Nov 03 15:30-1997 Dec 23 22:57 

P30805 1998 Jun 24 07:44-1998 Jun 26 09:05 

P40176 1999 Jan 04 06:25-2000 Feb 24 09:24 

P40177 1999 Jan 19 18:32-1999 Feb 01 20:15 



28 
132 

3 
57 

4 



116.6 
23.7 

133.7 
33.2 

473.3 
64.6 



S1/S2/GX1/GX2 
S1/S2/GX1/GX2 
S1/S2/GX1/GX2 
S1/S2/GX1/GX2 
S1/S2/GX1/GX2 
S1/S2/GX1/GX2 



total 



230 



845.1 



Number of observations conducted for each proposals. 

SI; Standard 1, S2; Standard 2, GX1; Good Xe-1 16s, GX2; Good Xe-2 16s. 

4 observations with large offset angle (Js 0.5 deg) were not used for the analysis. 
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Figure 1. Upper; long-term light curve of 3C 273. The count rates of the top layer from RXTE PCU 0/2 detectors are summed in the 
2.5—20 keV energy range. Lower; the count rates of All Sky Monitor (ASM) on-board RXTE. The data are binned in 10 day intervals. 



cause of the ASM problems described above). We have used 
a bin size of 5760 sec (approximately the RXTE orbital pe- 
riod) for the PCA plot. The 2PCU (PCUO/2) count rates 
are plotted for all data from 1996—2000. Expanded plots of 
the segments (a)-(d), indicated by arrows in Figure 1, are 
shown in Figure 2. An intensive hard X-ray monitoring for 
more than a month was conducted in 1996 simultaneously 
with IR telescopes (segment (b) in the Figures 1 and 2). Re- 
sults from this campaign have been discussed by McHardy 
et al. (1999). 

3 TEMPORAL STUDY 

3.1 Power Spectrum Density 

A power spectrum density (PSD) analysis is most commonly 
used to characterize the source variability. RXTE data pre- 
sented in this paper are the highest quality data ever re- 
ported in the X-ray band, enabling us to determine the PSD 
over a wider frequency range than attempted previously. An 



important issue, however, is the data gaps, which are un- 
avoidable for low-orbit X-ray satellites. Moreover, we must 
take great care since the observations were sparsely sched- 
uled; a net exposure of 845 ksec corresponds to only 0.7% 
of the total span of 4 years. In order to reduce the effect 
of data sampling, we follow a technique for calculating the 
PSD of unevenly sampled light curves (e.g., Hayashida et al. 
1998). 

The normalized PSD (NPSD) at frequency / is defined 

as 

P(/) = [" 2 (/) + fe 2 (/)-^ 2 tat Hr 

71-1 

a(f) = - VF J cos(2^ j ), 
n — ' 

j=0 
n-l 

Kf) = - V^sin^/tj), 
n — ' 

3=0 

(1) 
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Figure 2. Expanded plots of Figure 1. The light curves are sep- 
arately given for (a)-(d). The scale of time axis is same in all 
plots. 



where Fj is the source count rate at time tj (O^j^n — 1), 
T is the data length of the time series, and F av is the mean 
value of the source counting rate. The power due to photon- 
counting statistics is given by a^t ■ 

To calculate the NPSD of our data sets, we made light 
curves of three different bin sizes for RXTE data (256, 5760, 
and 43200 sec). We then divided each light curve into "seg- 
ments", which are defined as continuous parts of the light 
curve. If the light curve contains a time-gap larger than twice 
the bin size, we cut the light curve into two segments, one 
each side of the gap. We then calculate the power at fre- 
quencies / = k/T (1 ^ k ^ n/2) for each segment and 
take their average. The lowest frequency end (~ 10~ 6 Hz) is 
about half the inverse of the longest continuous segments, 
which is determined by the data in Figure 2(b) in our case. 

Figure 3 shows the NPSD calculated using this proce- 
dure. These NPSD data are binned in logarithmic intervals 
of ~0.2 (i.e., factors of approximately 1.6, as dictated by 
the data) to reduce the noise. The error bars represent the 
standard deviation of the average power in each rebinned 
frequency interval. The expected noise power due to count- 
ing statistics, <7g tat T /nF^ v , is shown in Figure 3 as a dashed 
line. As the NPSD has very steep power-law slope, with the 
NPSD decreasing as frequency increases, no detectable vari- 
ability exists in the light curve for / ^ 10~ 4 Hz (i.e., consis- 
tent with zero power for short time variability ^ 10 4 sec). 

A single power-law function is not a good representation 
of the NPSD, as the power-law fit above 10 -5 Hz is too steep 
for the data below 10 _s Hz (x 2 =14.7 for 8 dof; P( X 2 ) = 6%). 
A better fit was obtained using a broken power-law function, 
where the spectrum is harder below the break. The fitting 
function used was P(f) cc f~" L for / ^ f hr and P(f) cc 
f~ a for / ^ /br- The goodness of the fit was significantly 
improved; X 2 = 7.3 for 6 dof, P( X 2 ) = 30%. The best fit 
parameters are ql = 1.4 ± 0.2, a — 2.6 ± 0.1, and / br = 
(4.7 ± 1.5) x 10~ 6 Hz, respectively. This corresponds to a 
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Figure 3. Normalized PSD (NPSD) calculated from the light 
curves in Figure 1. The dotted line shows the best-fit broken 
power-law function. Full details are given in the text. 



variability time-scale for individual flares of f var ~ 3 days, 
as can be seen in Figure 2. 



3.2 Structure Function 

The structure function (hereafter SF) is a numerical tech- 
nique similar to the traditional PSD, but which has some ad- 
vantages for dealing with highly undersampled data. Since 
the SF is less affected by data gaps in light curves (e.g., 
Hughes, Alter, & Aller 1992), it may be a useful estimator 
for our study. The definitions of SFs and their properties are 
given in Simonetti, Cordes, & Heeschen (1985). The first or- 
der SF is defined as 



(2) 



where a(t) is a point in the time series (light curve) a, and 
the summation is made over all pairs separated in time by 
t. The term N is the number of such pairs. 

The SF is closely related to the PSD. If the structure 
function has a power- law form, SF{t) oc t p (f3 > 0), then the 
power spectrum has a distribution P(f) oc f~ a , where / is 
the frequency and a ~ (3+1. Obviously such an approximate 
relation is violated when (3 is close to zero, since both a and 
/3 approach zero for white noise. Nevertheless, the SF gives 
a crude but convenient estimate of the corresponding PSD 
distribution especially for red-noise type PSDs (e.g., Paltani 
et al. 1997; Cagnoni, Papadakis, & Fruscione 2001). 

In Figure 4, the SF is calculated from the light curve 
presented in Figure 1. We used all the data as a "continuous" 
observation in order to probe the variability on timescales 
as long as possible. RXTE (PCA) light curves binned in 
5760 sec intervals are used for the calculation. The resulting 
SFs are normalized by the square of the mean fluxes and 
are binned at logarithmically equal intervals. The measure- 
ment noise (Poisson errors associated with flux uncertainty) 
is subtracted as twice the square of Poisson errors on the 
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fluxes (and the noise level is given as the dashed line in the 
figure). 

As suggested by the NPSD in Figure 3, the resulting SFs 
are characterized by a steep increase (f3 ~ 1.5) in the time 
region of 0.1 ^ r ^ 3 days, above which the SFs show a sig- 
nificant roll-over (/? ~ 0.3). A gradual rise continues up to a 
time lag of r ~ 30 days, where a significant "bump" appears. 
This bump, and the "wiggling" features at the longest time- 
scales may be artefacts caused by the sparse sampling of 
the light curve (see Cagnoni, Papadakis, & Fruscione 2001; 
Kataoka et al. 2001 for a detailed discussion). The number 
of pairs (N in equation (2)) decreases as the time separa- 
tion t increases and, accordingly, the uncertainty becomes 
larger. 

The most rigorous study of the nature of this "bump" 
would be obtained by simulating the light-curves charac- 
terized by a certain SF, and filtering them with the same 
window as the actual observation. The resulting SF could 
then be compared with that adopted for the simulations. 
However, such an estimate is only possible when we already 
know the true "bump" structures of the system. The study 
of such structures is of interest in the light of claims for quasi 
periodic oscillation (QPO) on time-scales of ~ 100 day for 
other blazars (e.g., Rieger & Mannheim 2000; Abraham 
2000), but is beyond the scope of this paper. Here we limit 
discussion of the SF to time-scales shorter than 10 days, 
where N (see equation (2)) is large and the PSD and SF 
show good agreement. 

In Figure 4, we also show the SF of Mrk 421 calculated 
from ASCA observations over 5 years for comparison (solid 
line; reproduced from Kataoka et al. 2001). Mrk 421 is an 
"orthodox", and best studied blazar in the X-ray band. Two 
important characteristics are seen in the figure; (i) the struc- 
ture functions of both Mrk 421 and 3C 273 show a roll-over 
at 1 < t var < 10 days, and (ii) the normalized SF, which 
is proportional to the square of variability amplitude, is 1- 
2 orders of magnitude smaller for 3C 273 than for that of 
Mrk 421. 



4 SPECTRAL STUDY 

4.1 X-ray Spectral Evolution: power-law 
component 

Time variability of blazars is generally accompanied by sig- 
nificant spectral changes, which provide direct information 
about the acceleration, cooling and injection of electrons in 
the relativistic jet. Notably, such spectral information gives 
independent, and/or alternative information from the vari- 
ability studies presented in § 3. In this section, we investigate 
the photon spectra evolution of 3C 273 for each observation 
conducted from 1996 to 2000. We first divided the total ex- 
posure into one-orbit (5760 sec) intervals to investigate the 
most rapid evolution in photon spectra as possible. Spec- 
tral fits have been performed individually for 261 segments 
over the four years of observations. A single power-law func- 
tion and a photoelectric absorbing column Nn fixed at the 
Galactic value (1.79 x 10 20 cm~ 2 ; Dickey & Lockman 1990) 
represent most of the spectra well in the 2.5—20 keV data. 
The reduced x f° r 251 segments (out of 261) ranges from 
0.4 to 1.4 for 45 degrees of freedom, which corresponds to 
P(X 2 ) > 5%. 




ict 1 io° io 1 io 2 

time lag t (day) 



Figure 4. The structure function calculated from the light curves 
of Figure 1. The solid line shows the SF of the HBL Mrk 421 to 
compare with 3C 273 (QHB). Full details are given in the text. 

We summarize the results of power-law fitting in Fig- 
ure 5. The upper panel shows the changes in the 2—10 keV 
flux and the evolution of differential photon index (r) deter- 
mined in the 2.5—20 keV band. We are aware that RXTE 
is not sensitive to photons in the 2.0—2.5 keV. Neverthe- 
less, the 2—10 keV flux is specifically selected because most 
previous works have defined the X-ray flux in this energy 
band (e.g., Turner et al. 1990; Yaqoob et al. 1994; Haardt et 
al. 1998). The 2—10 keV flux changes dramatically during 
the observations in 1996-2000; from 2.25 x 10" 10 erg/cm 2 /s 
at the brightest (MJD 50400.8) to 0.62 x 10" 10 erg/cm 2 /s 
at the faintest (MJD 50621.7). During the 1996-1997 sea- 
son, the hard X-ray spectral photon indices stayed almost 
constant, with a differential photon index of T = 1.6 ± 0.1. 
Larger spectral variation is seen in the 1999—2000 season, 
where F varies from 1.5 to 1.9. 

The spectral evolution is more clearly seen in Fig- 
ure 6 where the correlation between the fluxes and the pho- 
ton indices are separately plotted for the 1996—1997 and 
1999-2000 seasons. For the 1996-1997 observations, the 
spectral shapes do not change significantly (r~1.6) despite 
large flux variations. The correlation factor between the 
flux and photon index is Rff = —0.19. A stronger correla- 
tion has been found for 1999—2000 observations, where the 
spectral index becomes softer when the source is brighter 
(Rft = -0.73). 

4.2 The Soft Excess and Broad Line Feature 

As mentioned above, 10 out of the 261 photon spectra can- 
not be fitted by a simple power-law function, in the sense 
that P(x 2 ) ^ 5%. There are two different types of departure 
from a power-law spectra; (i) a significant soft excess below 
4 keV, or (ii) a broad-line feature seen at 5—6 keV (in the 
observer's frame). 

Figure 7 shows examples of vF v photon spectra for each 
case. The left panel shows a concave feature with an appar- 
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Figure 5. The spectral evolutions of 3C 273 as a function of time, upper; changes in the 2—10 keV flux, lower; changes in the differential 
photon index V. 
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ent break at 4 keV (case (i)). A simple power-law fitting 
to the data gives \ 2 = 1-4 for 44 dof, where P(x 2 ) = 4%. 
Spectral fits with a broken power-law model significantly im- 
prove the goodness of the fit, as summarized in Table 2. Note 
that, below this break energy, the photon spectrum shows 
a very steep power-law index of Tl = 2.27 ± 0.19. A simi- 
lar soft excess has been found in the photon spectra of the 
segment closest in time (MJD 51246.15-51246.16), however, 
the statistical significance is relatively low and the simple 
power-law function gives a "good" fit {P(x 2 ) ^ 27%). 

The broad- line features, as illustrated in Figure 7 
(right), are observed 9 times during the 1996—1998 observa- 



tions. The results of spectral fits to a power-law plus Gaus- 
sian line model are given in Table 3. The line center energy 
is distributed between 5.8 and 7.1 keV in the quasar's rest 
frame. We cannot determine the exact shape and width of 
this broad line due to the limited sensitivity of the PCA 
on-board RXTE. Following Yaqoob & Serlemitsos (2000), 
we thus fixed the intrinsic line width at of c = 0.8 keV. The 
equivalent width of the line ranges from 108 to 210 eV al- 
though uncertainties are large. Taking these uncertainties 
into account, it seems premature to identify this structure 
as a broad Fe K a line only using the RXTE data. Compar- 
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Table 2. Spectral fits to 3C 273 with a broken power-law function. 



Obs.date 








2-10 keV Flux 


reduced \ 2 (dof) 


A X 2 |j 


(MJD) 




(keV) 




(10~ 10 erg/cm 2 /s) 




51246.10-12 


2.27±0.19 


4.02±0.37 


1.79±0.02 


1.17 


0.86(42) 


25.0 



a Photon index of the low-energy power-law component. 
^ The break energy. 

c Photon index of the high-energy power-law component. 

^ The decrease in \ 2 when the break is added to the simple power-law function. 



Table 3. Spectral fits to 3C 273 with a power-law function plus Gaussian line model. 



Obs.date 








EW^| 


2-10 keV Flux X 2 (dof) A X 2 ^ 


(MJD) 


(keV) 


(10- 4 /cm 2 /s) 


(eV) 


(10~ 10 erg/cm 2 /s) 



50281.14- 
50282.01- 
50282.14- 
50282.28- 
50282.35- 
50300.90- 
50400.82- 
50454.33- 
50988.32- 



•17 
05 
•17 
31 
38 
91 
83 
34 
36 



1.63±0.01 
1.62±0.01 
1.65±0.02 
1.62±0.01 
1.65±0.02 
1.62±0.02 
1.65±0.01 
1.63±0.02 
1.65±0.01 



6.32±0.19 
5.82±0.17 
6.10±0.30 
6.24±0.20 
6.43±0.19 
7.11±0.27 
6.05±0.23 
6.08±0.26 
6.10±0.21 



2.40±0.88 
2.91±0.98 
1.65±1.02 
2.44±0.90 
2.87±1.02 
2.53±1.23 
4.98±2.17 
5.44±2.65 
2.54±1.04 



187±68 

184±62 

108±67 

167±62 

206±73 

210±102 

132±58 

189±92 

117±48 



0.82 
0.88 
0.91 
0.92 
0.92 
0.93 
2.24 
1.74 
1.31 



0.75(43) 

0.78(43) 

1.18(42) 

0.85(43) 

1.21(43) 

1.02(300) 

1.04(43) 

1.09(43) 

1.09(43) 



30.8 
37.5 
11.7 
30.3 
32.1 
16.5 
21.4 
16.8 
25.5 



a Photon index of a power-law component. 

^ Line center energy. Intrinsic line width is fixed at 0.8 keV. 

c Intensity of a line. 

^ Equivalent width of a line. 

e The decrease in \ 2 when the Gaussian emission line is added to a simple power-law function. 
J Fitting was performed over 2.5-15 keV due to low photon statistics. 



ison with previous works, especially taken by instruments 
on-board other satellites, is necessary for confirmation. 

The broad line feature was firstly reported by Turner 
et al. (1990) using the data taken by EXOSAT and Ginga. 
RXTE detections of the Fe K a line from 3C 273 have been 
reported by Yaqoob & Serlemitsos (2000) for the observa- 
tions in July 17-18, 1996 (MJD 50181-50182) and June 
24-26, 1998 (MJD 50988-50990). Importantly, these obser- 
vations were conducted simultaneously with ASC'A. ASCA 
carried two different types of spectrometers, the Solid-state 
Imaging Spectrometer (SIS; Yamashita et al. 1997) and the 
Gas Imaging Spectrometer (GIS; Ohashi et al. 1996), both 
of which have much better line sensitivity than the RXTE 
PCA below 10 keV. Their results clearly indicate the exis- 
tence of a broad-line feature around 6 keV of EW ~ 100 eV 
(Figure 2 in Yaqoob & Serlemitsos 2000). Our analysis thus 
confirms the existence of broad line features for two obser- 
vations reported in Yaqoob & Serlemitsos (2000), and also 
adds new evidence from other observations in 1996—1997. 

It would be of interest to search for correlations be- 
tween the presence of the line and the 3C273 flux, however, 
the data available do not allow us to perform a convincing 
analysis of the line variability or of its relationship with the 
continuum variations. A future study along these lines will 
be very important for understanding where the cold matter 
emitting the fluorescence line is located with respect to the 
primary source. 



+ 



p 0.5 




3 5 7 10 20 30 3 5 7 10 20 30 

Photon Energy (keV) 

Figure 7. Examples of spectral fits to 3C 273 photon spectra 
when (a) a soft excess and (b) a broad line feature are present. 
The spectra are from data segments MJD 51246.10-51246.12 and 
MJD 50282.14-50282.17, respectively. 



4.3 Multi-frequency Spectrum 



In order to understand the origin of hard X-ray emis- 
sion from 3C 273 in more detail, we constructed a vL v 
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Table 4. Quadratic fit results of 3C 273 SED. 
components u p (Hz) L p (erg/s) L tot (crg/s) 

le IF 33 io 41 ^ iF 5 ^ 

heQ 10 20 io 46 - 2 10 47 1 

BBB IO 155 IO 462 IO 468 



a X-ray data taken during the brightest state are used for the 
fitting (MJD 50400.8; crosses in Figure 8.) 

mult i- frequency spectrum adding our new RXTE results. 
The resultant multi-frequency spectrum in the quasar's rest 
frame (z = 0.158) is shown in Figure 8. Crosses show 
RXTE data when the source was in the brightest state 
(MJD 50400.82-50400.83), whereas the filled circles show 
the data for the faintest state (MJD 50621.67-50621.69). 
Open circles and squares are archival data from literature; 
NED data base (radio to UV), ROSAT (soft X-ray; Staubert 
1992), OSSE, COMPTEL and EGRET on-board CGRO (7- 
ray; von Montigny et al. 1993; 1997). Three different peaks 
are apparent in the spectra; at low-energies (LE: radio to 
optical), high-energies (HE: X-ray to 7-ray), with the BBB 
(optical to soft X-ray) in between. 

To estimate the peak frequencies and luminosities of 
each component quantitatively, we fitted the spectral en- 
ergy distribution (SED) with separate polynomial functions 
of the form; \og{vL v ) — a + b logu + c (log!/) 2 , where a, 
b, c are constants. Similar polynominal fits were applied in 
Comastri, Molendi & Ghisellini (1995) for various blazars as- 
suming cubic functions. In the case of 3C 273, however, there 
are only negligible differences between the quadratic and 
cubic functions and hence, we adopt the simpler quadratic 
form in the fitting. The results of quadratic fits are sum- 
marized in Table 4, where L p is the peak luminosity in the 
vL v space, and Ltot is the integrated luminosity over that 
frequency range. 

It must be noted that the most of the data in Figure 8 
were not obtained simultaneously. It is clear from a number 
of observations conducted at various times and at various 
wavelengths that 3C 273 is variable from radio to 7-ray en- 
ergies (see von Montigny et al. 1997; Courvoisier 1998). Even 
for this particularly well-sampled and bright object, it has 
not yet been possible to obtain an actual "snap-shots" cov- 
ering the total energy band pass. Fortunately, the 30-year 
monitoring of 3C 273 reveals that the variability amplitude 
is factor of 4 at most energy bands (Tiirler et al. 1999). 
This is relatively small when plotted in the log(i/)-log(i/L„) 
plane, and we thus believe our discussion below is not af- 
fected significantly by the non-simultaneity of the data. We 
will comment on the effect of changes in physical parameters 
in §5.3. 



5 DISCUSSION 
5.1 Time Variability 

In §3, we have estimated the X-ray variability time-scale of 
3C 273 using two different methods. We found that a roll- 
over appears at a time-scale of t var ~3 day, which is similar 
to the X-ray variability time-scale of Mrk 421 (i V ar~ 1 day; 



Figure 4). However, a clear difference was found in the 
amplitude of the variability (normalization in Figure 4). In 
fact, the X-ray flux of Mrk 421 sometimes changes by more 
than a factor of 2 in a day (e.g., Takahashi et al. 2000), 
whereas flux variations of only 20—30% are detected for 
3C 273 (Figure 2). Two possibilities may be considered to 
account for such difference in variabilities. 

We first assume 3C 273 as a "pure" blazar, i.e., most of 
X-rays originate in beamed jet emission. Such an assump- 
tion may be valid since most of the time the X-ray photon 
spectra are well represented by a power-law function similar 
to other blazar-type objects. As many authors have sug- 
gested, variability in blazars may occur in a variety of ways, 
e.g., changes in acceleration rate of electrons, magnetic field 
strength and/or beaming factor. Most convincing scenario 
which may reproduce the observation well is that electrons 
are accelerated to higher energy during the flare than usual 
(e.g., Mastichiadis & Kirk 1997; Kataoka et al. 2000). If this 
is the case, the different variability amplitudes may result 
from different physical processes for the X-ray production. 

As expected from the multi-frequency spectra, the X- 
ray photons are thought to be emitted from the highest 
energy electrons via the synchrotron process for Mrk 421 
(Takahashi et al. 2000), whereas the inverse Compton emis- 
sion from the low-energy electrons are responsible for the 
case of 3C 273 (see Figure 8). Since the higher energy elec- 
trons cool faster (i CO oi oc 7 -1 , where 7 is the Lorentz factor 
of electrons; Rybicki & Lightman 1979), it is a natural con- 
sequence that both the temporal and spectral evolutions are 
most pronounced in the highest end of electron population, 
as observed in Mrk 421. This scenario may well explain the 
apparent difference of variability amplitudes in Mrk 421 and 
3C 273. 

Alternatively, we next assume 3C 273 as "half- 
Seyfert/half-blazar". It is clear that at some epochs at least, 
the X-ray flux has components that are not related to the 
beamed jet emission and that are observed to have fluxes 
similar to those discussed. This is particularly true of the 
soft excess and the broad line feature discussed in §4.2. Thus 
we may have observed a combination of non-thermal X-ray 
emission from the beamed jet, as well as the isotropic ther- 
mal emission from the accretion disk at the same time. The 
two types of emissions cannot be discriminated only from the 
X-ray data, however, variability in the lower-energy band, 
in particular the BBB, provides an important hint. 

The variability of the soft excess has been monitored 
for a long time. For example, Paltani, Courvoisier, & Wal- 
ter (1998) reported both optical and ultraviolet observations 
of 3C 273 covering the lifetime of the IUE satellite. They 
investigate the variability time-scales of BBB using a struc- 
ture function analysis. While resulting SFs are characterized 
by a steep increase as for the X-ray band (/3~1.5; Figure 4), 
a plateu appeared at much longer time scale (r ~ 0.5 year; 
but see Courvoisier et al. 1988 for an exceptionally rapid 
flare in the optical band). 

If t ~ 0.5 year is a "typical" variability time-scale 
of the non-beamed isotropical component in 3C 273, this 
would produce a "quasi-steady" underlying component to 
the "rapidly variable" beamed jet emission. Such an offset 
may significantly reduce the amplitude of X-ray variability 
originated in the jet, but the observed time-scales (i.e., du- 
ration of a flare) would not be changed (t~ 3 day). From 



RXTE Observations of 3C 273 between 1996 and 2000 '.Variability Timescale and Jet Power 9 

- , | , , , , | , , , , | , , , , - 
. 3C273 




10 15 20 

frequency log v (Hz) 



Figure 8. The multi-frequency spectra of 3C 273. Open circles and squares; archival data from the literature. Crosses; RXTE data 
when the source was the brightest, MJD 50400.8 (this work), filled circles; RXTE data when the source was the faintest, MJD 50621.7 
(this work). The solid line shows the SSC model for parameters B = 0.4 G, 5 = 10, 7max = 2000, and R = 2.5 X 10 16 cm. The dashed 
line shows a model for parameters B = 0.4 G, 8 = 10, 7max = 2000, and R = 5.0 X 10 16 cm with moderate electron number density. All 
parameters are determined self-consistcntly as discussed in the text. The dotted line shows a quadratic fit to the BBB to guide the eye. 



this view-point, the difference of the variability amplitude 
in Mrk 421 and 3C 273 may be explained by the impor- 
tance of underlying isotropic component over the beamed 
jet emission. 

5.2 Spectral Evolution 

In §4, we showed that most 3C 273 photon spectra ob- 
served between 1996 and 2000 can be represented by a sim- 
ple power-law function with T = 1.6 ± 0.1. This is consis- 
tent with previously published results. Turner et al. (1990) 
suggest that the spectral index of 3C 273 in the 2—10 keV 
band is slightly harder than the "canonical" AGN index of 
T = 1.7, though the difference seems to be very small (1.5 
compared to 1.7). Such a difference, if confirmed, may be 
due to a different physical origin of X-ray production. 

In fact, X-ray photons are thought to be emitted from 
the Compton scattering of disk photons for Seyfert-type 
AGNs, whereas accelerated electrons in the jet produce non- 
thermal emission in the case of blazars. In the standard 
shock scenario (Fermi acceleration), electrons have a pop- 
ulation AT (7) cc 7~ 2 , where AT (7) is the number density of 
electrons per energy. Such accelerated electrons emit pho- 
tons with a spectrum of a form cx v~ , through both the 
synchrotron and inverse-Compton processes (e.g., Rybicki 
& Lightman 1979). This is very close to the observational 
results. 

If the hard X-ray photon index reflects the power-law 
index of the accelerated electrons (s of 7 _s ), only small spec- 



tral changes would be expected in the photon spectra cor- 
responding to the low-energy end of electron populations 
(as for the X-ray photons of 3C 273). This is because low 
energy electrons do not cool before leaving the emitting re- 
gion, and hence keep their "original" (injected) spectral in- 
formation. As suggested by Kirk, Rieger, & Mastichiadis 
(1998), changes in the electron injection rate could produce 
flux variations, but the photon spectra do not evolve signif- 

-(s + l) 

icantly (oc v 2 ). Such a scenario may well explain the 
spectral evolution observed in the 1996—1997 season, where 
the photon index stays almost constant (T~1.6) regardless 
of the source flux variations. 

In 1999—2000, a different spectral evolution was found, 
with the spectra becoming steeper (up to T~ 1.8) as 
the source becomes brighter. This may be related to the 
underlying isotropic component emitted from the accretion 
disk or corona, as we have discussed in § 5.1. It is interesting 
to note that the exceptional "concave" photon spectrum was 
found only in this season (Figure 7 left). Such excess emis- 
sion may be the hard-tail of the BBB, as suggested from a 
ROSAT observation. Staubert (1992) reported a very steep 
X-ray spectrum of Tl ~ 2.6 in the ROSAT bandpass (see 
also Figure 8), which is marginally consistent with our re- 
sults (Fl = 2.3 ±0.2; Table 2). It is still unclear if the excess 
emission is really the hard tail of the BBB, as this feature 
was only clearly observed once, and RXTE has relatively 
poor sensitivity in the soft X-ray band. However, it seems 
plausible that non-beamed, thermal emission and its repro- 
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cessed flux becomes more important in 1999-2000 season, 
compared to those in 1996-1998: future studies are neces- 
sary to confirm this. 

The broad line feature is probably the Fe fluorescent 
line, as already discussed in Turner et al. (1990; EXOSAT 
and Ginga), Cappi et al. (1998; ASCA) and Yaqoob & Ser- 
lemitsos (2000; ASCA and RXTE). Interestingly, the Fe 
line is usually not visible in the RXTE observations; only 
9 of 261 observations clearly show line features, which is 
at odds with any thought that the hard X-ray photons are 
emitted by the same process with those in Seyfert galax- 
ies. This suggests again that the X-ray photons of 3C 273 
are ordinarily emitted from the relativistic jet, but that 
occasionally the emission from the accretion disk is visible. 
(Although it should be kept in mind that the RXTE PCA 
has a relatively low sensitivity for the broad line feature, 
and that observations with a more sensitive spectrometer 
may have detected the line more frequently.) The visibility 
of the Fe fluorescent line may depend on the balance between 
the jet emission, blue bump, and/or Compton scattering of 
disk photons. Future studies along these lines may be able to 
discriminate between various models of disk/jet connection. 

In summary, we conclude that the hard X-ray emission 
of 3C 273 may occasionally contain a significant contribution 
from the accretion disk, but most hard X-rays are likely to 
originate in inverse Compton radiation from the relativistic 
jet. 

5.3 Multi-frequency spectrum 

Multi-frequency spectra provide information of physical 
quantities relevant for jet emission, e.g., the magnetic field, 
the size of the emission region, the maximum energy and the 
density of relativistic electrons. Various authors have deter- 
mined jet parameters in a self-consistent manner for blazars 
(e.g., Mastichiadis & Kirk 1997; Kataoka et al. 2000; Kino, 
Takahara, & Kusunose 2002), but only a few attempts have 
been made for 3C 273. We note again that part of the emis- 
sion from 3C 273 may not be related to the beamed emission 
as discussed above. Nevertheless, we will consider this ob- 
ject as a "blazar" for the following reasons; 

(1) Contrary to many Seyfert galaxies (e.g., Mushotzky, 
Done & Pounds 1993), there is no hint of a "reflection hump" 
corresponding to the Compton reflection of primary compo- 
nent. 

(2) Although some exceptions have been found, the X-ray 
photon spectra are generally well represented by a power- 
law function which smoothly connect to the 7-ray bands, as 
for other blazar-type objects. 

Von Montigny et al. (1997) examined three models to 
reproduce the multi-frequency spectrum of 3C 273: (i) the 
synchrotron self-Compton model (SSC model; e.g., Inoue & 
Takahara 1996), (ii) the external radiation Compton model 
(ERC model; e.g., Sikora, Begelman, & Rees 1994), and 
(iii) the proton induced cascade model (PIC model; e.g., 
Mannheim & Biermann 1992). Although 3C 273 is the one 
of the best sampled and best studied objects across the en- 
tire electromagnetic spectrum, they conclude that the data 
were still insufficient to discriminate between these models. 

In the following, we thus adopt the simple one-zone 
SSC model to describe the spectrum. The resultant phys- 
ical quantities provide a lot of feed-back to test the validity 



of our assumptions. We will also comment on the applicabil- 
ity of other emission models, especially the ERC model, in 
the appendix. We do not consider the PIC model, in which 
the jet is composed of e~-p plasma (baryonic) rather than 
e~-e + plasma (leptonic). Such a situation might be possible, 
but observationally, there are no reasons for considering ad- 
ditional, baryonic emission to account for the overall spectra 
(see Sikora & Madejski 2001 for a detailed discussion). 

The jet parameters are tightly connected with the ob- 
served quantities, particularly the v p and Ltot of each emis- 
sion component (see §4.3). In the simple model we adopt 
here, the radiation is due to a homogeneous jet component 
moving with a bulk Lorentz factor Tblk = (1 — ft 2 )^ 1 ^ 2 at 
an angle to the line of sight of 9 ~ 1/Fblk- 

We first assume that the peak emission (L p ) of the low- 
energy (LE) synchrotron component and the high-energy 
(HE) inverse-Compton (SSC) components arise from the 
same electron population with a Lorentz factor 7 P . The peak 
frequencies of LE and HE are related by 

^HE,p = g7p^LE,p, (3) 

where fLE, P = 10 13 ' 5 Hz and i^je.p = 10 20 Hz, respectively 
(Table 4). We thus obtain 7 P = 2.0xl0 3 . The synchrotron 
peak frequency, i^le, p , is given by 

^ LE , P ~3.7x 10 6 S 7 2 ^- Hz, (4) 

where B is the magnetic field strength and S — rg LK (l — ft 
cos #) _1 ~ Tblk is the relativistic beaming (Doppler) factor. 
The magnetic field B and its field densities are thus derived 

B = 0.4 x (-^) G, U B = 6.7 x 10~ 3 (-^) 2 erg/cm 3 . (5) 


Note that VLBI observations of 3C273 set a lower limit for 
the Lorentz factor of Tblk > 10, implying that 5 ~ Tblk 
~ 10 if the angle to the line of sight of #~ 1/Fblk (e.g., 
Pearson et al. 1981; Krichbaum et al. 2002). 

The ratio of the synchrotron luminosity, Z/LE,tot, to the 
inverse Compton (SSC) luminosity, LHE.tot is 

£<HE,tot U Bync . . 

7 = It' (6) 

-tJLE.tot UB 

where U sy nc is the synchrotron photon energy density in 
the co-moving frame of the jet. Using the observed value of 
L HB ,tot = 10 47 ' 1 erg/s and L LE , tot = 10 46 ' 7 erg/s (Table 4), 
we obtain 

[/ sync = 2.0x 10~ 2 (-^) 2 erg/cm 3 . (7) 


Assuming a spherical geometry for the emission region, 
the synchrotron luminosity is given as 

iLE.tot = 4nR 2 c5 4 Usync erg/s, (8) 

where R is the radius of the emission region. This radius is 
estimated as 

i? = 2.5x 10 16 (^) cm. (9) 


Remarkably, the size of the region derived here corre- 
sponds to the variability time-scale of t var ~ R/(cS) ~ 1 day, 
qualitatively consistent with those expected from the tem- 
poral studies in §3. Also, if the jet is collimated within a 
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cone of constant opening angle 9 ~ 1/Fblk, and the line- 
of-sight extent of the shock is comparable with the angular 
extent of the jet, one expects that the X-ray emission site is 
located at distance D ~ io 17 ~ 18 cm from the base of the 
jet. This implies that the radiation discussed here is emitted 
from the inner-most part of the jet, i.e., "sub-pc-scale" jet. 
(This must be clearly distinguished from the "large scale" 
jet described later in this section.) 

The SSC spectrum is self -consistently calculated as- 
suming the parameters derived above. We assume an elec- 
tron population of the form AT (7) oc 7~ s exp( — 7/7 P ), where 
s is set to be 2. The solid line in Figure 8 shows the over-all 
model for the parameters R — 2.5 x 10 16 cm, B = 0.4 G, 
8 — 10, 7 P = 2.0 x 10 , while the dashed line shows the case 
when the region size is larger by a factor of 2 with moderate 
electron number density (and other parameters unchanged) . 
The whole spectrum is adequately represented, except for 
the discrepancy in the radio band. Such discrepancies are 
common for one- zone model fitting of blazar SEDs, as one- 
zone models cannot account for the low-energy emission, 
which is thought to be produced in a much larger region of 
the source (e.g., Marscher 1980). 

The electron energy density, U e , is calculated to be 6.3 x 
10 _1 erg/cm , where we set the minimum energy of electrons 
to be 7min = 1 • The kinetic power of the electrons which emit 
the observed photons is thus estimated to be 

L kin ~ 7ri? 2 C r| LK f/e ^ 3.7 xl0 45 (^) 2 erg/s. (10) 

It is interesting to compare these radiative/kinetic lu- 
minosities to discuss the "power balance" between the accre- 
tion disk and the jet (see also Celotti, Padovani, & Ghisellini 
1997; Celotti, Ghisellini, & Chiaberge 2001 for the theoreti- 
cal approach). Assuming that the BBB is emitted isotropi- 
cally from the accretion disk, the electro-magnetic luminos- 
ity in the accretion disk is Ldisk — £<BBB,tot = 10 46 ' 8 erg/s 
(Table 4). The "radiative power" and "Poynting power" of 
the jet in the observer's frame are given as 

L rad ~ TvR 2 cT 2 BLK U pb ~ 5.0 x 10 44 (-i^) 2 erg/s (11) 



L B ~7ri? 2 C r BLK ir B ^3.9xl0 43 (^) 2 erg/s, (12) 



respectively, where U pn is the photon energy density in the 
jet co-moving frame; U p h — U Byn c + f/ssc- By comparing 
idisk, irad, Lb, and Lkin, we derive the following properties 
of the "power balance" in sub-pc-scale jet; 

(i) the kinetic power of relativistic electrons amounts to 
only 5—10% of the BBB emission, which is thought to be 
radiated in the vicinity of central black-hole. 

(ii) the magnitudes of powers are Lkin ^ L r a,d ^ Lb- 

We finally comment on the large scale X-ray jet ob- 
served by Einstein (Harris & Stern 1987), ROSAT (Roser 
et al. 2000) and Chandra (Sambruna et al. 2001; Marshall et 
al. 2001). The X-ray jet is ~ 8" (1" = 2.4 kpc at the distance 
of 3C 273) long and has a knotty morphology, starting with 
a bright, resolved knot ~ 13" from the core. This "kpc-scale 
jet" is much fainter than the "sub-pc-scale jet" . Indeed, the 
integrated X-ray flux of the kpc-scale jet is 6.9±0.6 njy at 
1 keV (Marshall et al. 2001), which corresponds to ~0.1% 
of the X-ray flux emitted from the "sub-pc-scale jet" . The 
spectral energy distributions from radio to X-rays of the in- 
dividual knots show a variety of shapes. Sambruna et al. 



(2001) fit the multi-frequency data of each knot with the 
ERC/CMB model, where the cosmic microwave background 
(CMB) photons are up-scattered into the X-ray band. For 
example, they obtain the best fit parameters for region A 
(see Sambruna et al. 2001 for definition) of electron index 
s — 2.6, normalization K = 8.1 x 10~ 3 cm~ 3 , magnetic field 
B = 1.9 x 10" 6 G, region size R = 5 x 10 cm, and Doppler 
factor 8 = 5.2. 

Assuming these quantities, we estimate the ki- 
netic power in the 10 kpc scale jet as L k in,iok P c = 
1q47.i/_20\0.6 er g/ s which is about 2 orders of magni- 

7min 

tude larger than that in the sub-pc-scale jet discussed above 
(£/kin,sub- P c ~ 4xl0 45 erg/s; see equation (10)). It should be 
noted that the kinetic power of electrons strongly depends 
on both the minimum energy (7mm) and the power-law in- 
dex (s) of the assumed electron population. Part of the dif- 
ference may be explained by different parameters selected 
above, however, it seems difficult to completely explain the 
discrepancy by this fact alone. In fact, if we assume 7 m i n = 1 
for the 10 kpc scale jet, as was assumed for the sub-pc-scale 
jet, Lkin.iokpc becomes larger, increasing the discrepancy. 
Similarly, an electron population with s — 2.0 would mini- 
mize Lkin.iokpc, but such parameters would not fit the ob- 
served SED so well. Future observations of the 10 kpc scale 
jet with wider energy bands are eagerly awaited. 

If the difference between the sub-pc and the 10 kpc 
scale jet powers is confirmed, this may imply that ~100 
times the "visible" kinetic energy is hidden at the bottom of 
the jet. It should be noted that the kinetic power discussed 
above takes account of only relativistic electrons, in other 
words, the contribution from thermal electrons as well as 
protons (either relativistic or cold) have completely been 
neglected. These non-relativistic electrons and/or protons 
can not contribute to the emission in sub-pc-scale jet. 

Only a small fraction of them, probably less than 1 % 
in number, may be "picked up" by the shock acceleration 
process, to add to the relativistic electron population. Such 
a low acceleration efficiency may well be understood by the 
so-called internal shock model as discussed in Spada et al. 
(2001) and Kataoka et al. (2001). Further discussion along 
this line is now in progress. For example, Tanihata (2002) 
suggests that, in order to explain the observed variability 
properties, the velocity difference of the two colliding shells 
must be rather small. For this case, the efficiency is esti- 
mated to be as small as < 0.01 %. 

The enormous kinetic energy could be released for the 
first time at large distances (e.g., ^ 10 kpc) by completely 
different mechanisms of energy dissipation. One such pos- 
sibility is the external shock model, wherein shocks arise 
when outflowing jet plasma decelerates upon interaction 
with dense gas clouds originating outside the jet (e.g., Der- 
mer & Chiang 1998). The precise nature of the required 
gas clouds is uncertain, but shock acceleration takes place 
more efficiently than those for the internal shock model. 
Meanwhile, the different jet powers may imply that part 
of the kinetic power carried by protons in the sub-pc-scale 
jet is efficiently transferred to the electrons in 10 kpc scale 
jet, though the mechanisms of energy dissipation are com- 
pletely uncertain. Future deep observations of large scale 
jets by Chandra and Newton may be able to confirm our 
suggestions about the power consumption in the 3C 273 jet. 
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6 CONCLUSION 

We have analyzed the archival RXTE data available for 
3C 273 between 1996 and 2000. A total of 230 observa- 
tions amounts to the net exposure of 845 ksec over the 4 
years. This is the longest, and most densely sampled, expo- 
sure for this object in the hard X-ray band. Both the PSD 
and the SF show a roll-over with a time-scale on the or- 
der of ~ 3 days, although the lower- frequency (i.e., longer 
time-scale) variability is still unclear. We found that the 
variability time-scale of 3C 273 is similar to those observed 
in TeV gamma-ray emitting blazars, whereas the variabil- 
ity amplitude is an order of magnitude smaller. Considering 
that the hard X-ray spectra of 3C 273 generally maintains 
a constant power- law shape with T ~ 1.6 ± 0.1, beamed, 
inverse Compton emission inside the jet is the most likely 
the origin of the X-ray /7-ray emission. Two kinds of excep- 
tions have been found, which may be interpreted either as 
the hard-tail of the BBB or the emission from the accre- 
tion disk occasionally superposed on the jet emission. From 
a multi-frequency analysis, we constrain the physical quan- 
tities relevant for the jet emission. We argue that, (i) the 
kinetic power carried by relativistic electrons corresponds 
to only 5—10% of the disk luminosity, and (ii) the various 
powers in the sub-pc-scale jet are ranked Lkin ^ £rad ^ Lb- 
The connection between the sub-pc-scale jet and the 10 kpc 
scale jet remains uncertain, but our work suggests that the 
most of the jet power might be hidden at the base of the 
sub-pc-scale jet, and effectively released at the 10 kpc scale 
via a completely different mechanism of energy dissipation. 



APPENDIX: APPLICABILITY AND VALIDITY 
OF THE SSC/ERC MODELS 

Finally, we comment on an alternative scenario which may 
account for the overall spectra of 3C 273, namely the exter- 
nal radiation Compton (ERC) model. 

A number of quasars show complicated multi-frequency 
spectra which cannot be readily fitted with a simple SSC 
model. This is mostly because the 7-ray flux strongly dom- 
inates the radiative output, and the 7-ray spectra are well 
above the extrapolation of the X-ray spectra (e.g., Ghisellini 
et al. 1998). One explanation is that the SSC process domi- 
nates in the X-ray range, while the ERC process dominates 
in 7-rays (e.g., Inoue & Takahara 1996). Such a disconti- 
nuity/complexity cannot be seen for 3C 273 in Figure 8, 
suggesting that both the X-rays and 7-rays are produced by 
the same radiation process, i.e., the SSC process - that is 
the main reason of why we apply the SSC model to the data 
in § 5.3. 

However, this explanation is unable to reject other emis- 
sion models, in particular the ERC model. This is because 
most of data are obtained non-simultaneously, and we do not 
have an actual "snap-shot" of the overall emission. In fact, 
the ERC model provides similarly good fit for the case of 
3C 273, as discussed in detail in von Montigny et al. (1997). 
In this appendix, we consider in what situation the SSC 
dominates the ERC process, and what the expected impact 
would be if the ERC model is responsible for the overall 
emission from 3C 273. 

Within the framework of ERC models, a number of pos- 
sibilities exist for the origin of the seed photons which are 



up-scattered to 7-rays. Dermer & Schlickeiser (1993) sug- 
gest that the direct emission from the central core irradiates 
the emission region inside the jet, which may provide a suf- 
ficient number of seed photons. In the observer's frame, the 
isotropic luminosity of the central core is approximated by 



i disk ~ L B BB,tot = 6.3 x 10 erg/s 



(13) 



(see Table 4). In the co- moving frame of a jet with bulk 
Lorentz factor Fblk, the photon energy density which is 
produced by the core emission is given as 

L C orc F BLK (1 + P — 2/3) L COFC 1 



u, 



47ni 2 



4ttcP 4r| LK < 



(14) 



where d is the distance from the central core. By normalizing 
by typical values of Tblk and d, we obtain 



, „ i/O.l pc. 2/ 10 

4.6 x 10 3 ( — j—) 2 {- 



' Tblk ' 



erg/cm 3 



(15) 



Comparing this with the synchrotron photon energy density, 
[/ sync = 2.0xl0~ 2 (10/5) 2 erg/cm 3 (see equation (7)), direct 
emission from the central core cannot be a dominant source 
of seed photons for the ERC process. This is because the jet 
plasma is moving away from the central core at highly rela- 
tivistic speed (Fblk— 10), and hence is strongly red-shifted 
as measured in the source co-moving frame. 

The second candidate source of the external seed pho- 
tons is the "diffuse" radiation field, i.e., a radiation field 
having a significant non-radial component at large distance 
from the central engine. Such a radiation component may be 
produced by scattering or reprocessing of a portion of the 
central radiation by irradiated clouds and / or the inter-cloud 
medium (Sikora, Begelman, & Rees 1994). Importantly, all 
photons directed inward toward the central source are blue- 
shifted, and enhanced by a factor of Tblk in the frame co- 
moving with the jet. Assuming that a certain fraction r (< 
1) of the central luminosity is reprocessed into such a diffuse 
radiation field, we obtain 



U, 



tL, 



ext,out — 



core J- BLK 



47Td 2 



(16) 



If the SSC process is more likely origin of the X-ray/7- 
ray emissions, as we have assumed in § 5.3, inverse Comp- 
tonization of external photons does not dominate the syn- 
chrotron photons and U syn c ^ (Text, out must hold. This limits 
the fraction r very tightly to 



T 1.1 X 10~*( 



d 



0.1 pc ; Tblk 



(17) 



meaning that less than 0.01 % of the BBB luminosity could 
be reprocessed into the "isotropic" external radiation field 
in the case of 3C 273. One difficulty in understanding the 
overall photon spectrum of 3C 273 with the SSC model is 
why such a low efficiency is achieved in this object - unfor- 
tunately, we still do not have an answer to this question. 

One may therefore consider the opposite extreme, when 
the X-ray and 7-ray emissions are dominated by the ERC 
process, rather than SSC emission. In such a situation, the 
ratio of the synchrotron luminosity, Z/LE,tot, to the inverse 
Compton (ERC) luminosity, LHEtot, is (see also equation 
(6)) 



:,t0t out 



(18) 
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The magnetic field density is estimated to be 



Ul 



47rd 2 



(19) 



The reprocessing factor r is very uncertain, but it must be 
larger than 10" 4 if the ERC process dominates (see equation 
(17)). Normalizing r with 10~ 3 , and assuming the observed 
value of I/LE,tot, iHE.tot, and L CO r C — ^BBB,tot (Table 3), we 
obtain 



B ~ 1.4 x ( 



10 



-3> ( 1Q 



d 



(20) 



Thus in the framework of the ERC model, the overall emis- 
sion of 3C 273 may be reproduced if the magnetic field 
strength is a bit stronger than the SSC case. 

In both the SSC and ERC models, the low energy com- 
ponent is thought to be produced by the same, synchrotron 
process. The synchrotron luminosity, Z/LE,tot behaves as 



oc UbUzR* 



(21) 



In order to explain the observed LLE.tot with a different 
magnetic field strength B (~0.4 G for the SSC and ~1.4 G 
for the ERC model), the electron energy density U e must 
be varied since R is constrained to 10 16-17 cm from the ob- 
served variability time-scale. Thus the U e of the ERC model 
could be smaller by about an order of magnitude than we 
have assumed for the SSC model, U e ~ 5xl0~ 2 erg/cm 3 . 
This may lead to an over-estimate of the kinetic power if 
the overall X-ray/7-ray spectra are dominated by the ERC 
(see equation (10)). However, what is important is that as- 
suming either the SSC or the ERC model, our discussion 
and conclusions about the power balance in the jet (see § 
5.3) are not affected significantly. 
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